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Abstract 

The  paper  shows  and  discusses  a  procedure  of  parameter  estimation  applied  to  the  evaluation  of  some  operating  parameters  of  a  proton-exchange 
membrane  fuel  cell  (PEMFC).  First,  a  brief  literature  review  about  the  main  parameters  (exchange  current  density,  cell  resistance,  internal  current 
density  and  limiting  current  density)  has  been  done.  Then  the  analytical  model  adopted  to  describe  the  polarization  curve  has  been  discussed. 
Based  on  this  model,  a  parameter  analysis  has  been  done,  and  it  has  been  shown  that  three  parameters  of  the  cell  polarization  curve  model  can  be 
simultaneously  estimated:  the  cathode  exchange  current  density,  the  cell  resistance  and  the  internal  current  density.  To  evaluate  these  parameters  both 
a  set  of  our  measurements  on  a  PEM  single  cell  (active  area  of  25  cm2  and  Nafion  115  membrane)  and  data  from  other  authors  has  been  considered. 

The  cell  has  been  fed  with  pure  hydrogen  and  air,  the  cell  temperature  has  been  varied  from  50  °C  to  80  °C,  and  accordingly  the  reactants  have 
been  introduced  in  the  cell  humidified  at  the  same  temperature.  The  parameters  have  been  estimated  in  each  operating  conditions  of  the  cell,  and 
their  behavior,  as  a  function  of  the  cell  operating  temperature,  has  been  discussed. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction  and  aims  of  the  paper 

The  performance  of  a  proton-exchange  membrane  fuel  cell 
(PEMFC)  can  be  expressed  through  the  analytical  formulation 
of  the  polarization  curve  [1,2];  different  approaches  could  be 
considered,  and  the  models  formulation  introduces  some  oper¬ 
ating  parameters.  Among  the  most  important  parameters  are: 
exchange  current  density  at  the  anode  z'o,a  and  cathode  z'o,c  elec¬ 
trodes;  cell  resistance  r;  internal  current  density  /n;  limiting 
current  density  at  the  anode  i\tSL  and  cathode  i\?c  electrodes.  In 
literature,  the  discussion  about  the  simultaneous  estimation  of 
the  previous  parameters  is  not  particularly  developed.  Besides 
that,  there  are  few  evaluations  concerning  the  behavior  of  the 
cell  parameters  versus  cell  operating  temperature.  Therefore, 
this  paper  starts  with  an  analysis  about  the  available  literature 
concerning  the  cell  parameters  evaluation.  Then,  a  regression 
model  of  the  polarization  curve,  adopted  for  the  parameter  esti¬ 
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mation  procedure,  is  described,  and  the  characteristic  parameters 
are  outlined.  After,  the  parameter  estimation  procedure  (inverse 
model  [3])  is  described,  which  makes  use  of  the  experimental 
data  of  the  independent  variables  (such  as  current  density  and 
operating  temperature)  and  the  dependent  variables  (such  as  cell 
voltage).  Some  preliminary  tests  made  on  the  model  showed 
the  number  and  type  of  parameters  which  had  been  possible  to 
correctly  estimate  from  the  experimental  data.  After,  the  main 
parameters  of  the  polarization  curve  of  a  single  PEM  have  been 
estimated  (a  brief  description  of  the  experimental  setup  and  of 
the  fuel  cell  adopted  is  also  provided).  Then,  the  behavior  of 
the  estimated  parameters  as  a  function  of  an  independent  vari¬ 
able  (cell  operating  temperature)  has  been  discussed.  Finally, 
the  parameter  estimation  procedure  has  been  extended  to  other 
sets  of  experimental  data  and  the  agreements  and  discrepancies 
are  discussed. 

2.  Brief  literature  analysis 

A  wide  and  detailed  discussion  about  these  parameters  is 
available  in  literature. 
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Nomenclature 

a 

activity  of  a  substance 

Bp 

backpressure  of  reactants  (bar) 

c 

concentration  of  a  chemical  species  (molm-3) 

Coo 

concentration  of  a  chemical  species  before  the  dif¬ 
fusion  layer  (mol  m-3) 

Cl 

parameter  of  the  cell  resistance  expression 
(£2  cm2) 

C2 

parameter  of  the  cell  resistance  expression 
(£2  cm4  A-1) 

D 

diffusivity  of  a  chemical  species  (cm2  s-1) 

E 

open  circuit  voltage  of  the  single  PEM  fuel  cell 
(V) 

F 

Faraday  number  (96,485  C  mol-1) 

g 

molar  Gibbs  free  energy  (J  mol-1) 

i 

current  density  (A  cm-2) 

h 

limiting  current  density  in  a  PEM  fuel  cell 
(A  cm-2) 

h,  a 

anode  limiting  current  density  (A  cm-2) 

h,  c 

cathode  limiting  current  density  (A  cm-2) 

in 

internal  current  density 

fi),a 

anode  exchange  current  density  (A  cm-2) 

fi),c 

cathode  exchange  current  density  (A  cm-2) 

fi),el 

exchange  current  density  in  an  electrode 
(A  cm-2) 

/ 

current  (A) 

M 

molar  mass  (g  mol-1) 

n 

number  of  electrons  participating  in  a  reaction 

nd 

electro-osmotic  drag  coefficient  (moln2o  mol“+) 

N 

molar  flux  per  surface  area  (mol  s-1  cm-2) 

P 

partial  pressure  of  a  gas  (bar) 

PEMFC 

proton-exchange  membrane  fuel  cell 

r 

cell  resistance  of  the  PEM  fuel  cell  (Q  cm2) 

rel 

electronic  resistance  (£2  cm) 

^ion 

ionic  resistance  (£2  cm2) 

R 

universal  gas  constant  (J  mol-1  K-1) 

fin 

membrane  thickness  (cm) 

T 

temperature  (°C) 

rr S 

1  a 

temperature  of  saturated  inlet  anode  flow  (°C) 

rrS 

1  c 

temperature  of  saturated  inlet  cathode  flow  (°C) 

Vc 

voltage  of  the  single  PEM  fuel  cell  (V) 

Greek  letters 

aa 

anode  transfer  coefficient 

ac 

cathode  transfer  coefficient 

af 

transfer  coefficient  on  an  electrode 

a? 

transfer  coefficient  on  an  electrode 

8 

Nemst  diffusion  layer  thickness  (cm) 

Sjt 

percolation  threshold 

y 

empirical  constant  of  diffusivity  correction 

0 

overpotential  (V) 

0  act 

activation  overpotential  (V) 

7?act,a 

anode  activation  overpotential  (V) 

7?act,c 

cathode  activation  overpotential  (V) 

0  cone 

concentration  overpotential  (V) 

7?conc,a 

anode  concentration  overpotential  (V) 

0conc,c 

cathode  concentration  overpotential  (V) 

7?  ohm 

ohmic  overpotential  (V) 

X 

degree  of  humidification  of  the  membrane 
(molH2o  mol~Q3_) 

■^h2 

hydrogen  excess 

^air 

air  excess 

a 

membrane  ionic  conductivity  (£2  cm)-1 

2.1.  Exchange  current  density 


It  pertains  to  the  specific  electrochemical  reaction,  and  it 
is  mainly  a  function  of  the  electrode  characteristics:  type  and 
quantity  of  catalyst,  dimension  and  distribution  of  the  catalyst 
particles,  active  surface.  Considering  the  operation  variables,  it 
is  linked  especially  to  the  cell  temperature. 

The  catalyst  type  has  an  important  role.  Some  indicatives 
values,  related  to  the  anode  electrode  of  a  PEM  at  T-  25  °C: 
Ag,  /o,a  =  4  x  10-4  A  cm-2;  Ni,  /o,a  =  6  x  10-3  A  cm-2;  Pt, 
z'o,a  =  5  x  10"1  cm-2  [4]. 

The  increase  of  the  operation  temperature  would  cause  an 
increase  of  the  exchange  current  density,  because  a  higher  tem¬ 
perature  value  allows  an  improvement  of  the  reaction  activation 
[4]. 

In  a  PEM  fuel  cell,  the  value  of  the  exchange  current  density 
at  the  cathode  electrode  is  considerably  low  compared  to  the 
value  at  the  anode  electrode,  and  therefore  the  anodic  activa¬ 
tion  overvoltage  are  usually  negligible.  As  an  example,  typical 
values  are:  /o,a  =  0.2  A  cm-2  and  io,c  =  1  x  10-4  A  cm-2  [4].  In 
ref.  [5],  a  value  of  /o  =  1  x  10-4  A  cm-2  is  considered,  even  if 
it  is  not  explained  whether  at  the  anode  or  (probably)  at  the 
cathode  side.  In  ref.  [6],  the  values  are:  /o,a  =  0.0538  A  cm-2 
and  z’o,c  =  1.0764  x  10-6  A  cm-2  at  T=  25  °C;  in  the  paper,  an 
analytical  expression  of  the  exchange  current  density,  taken 
from  Berger  [7]  and  associated  to  a  generic  electrode  reaction 
a  A  +  nt~  £>B,  is  also  described: 


/0  =  n  •  F  •  k° 


exp 


—  A  Fq\ 
RTJ 


a-(\—  ael) 
A 


C 


b-a el 
CB 


where  A  FQ  is  defined  as  free  standard  energy  of  activation,  which 
is  finally  the  change  of  the  Gibbs  free  energy  of  reaction,  while  k° 
is  a  parameter  connected  to  the  reaction  speed.  In  ref.  [8],  there 
are  two  different  analytical  expressions  related  to  the  anode  and 
cathode  exchange  current  density: 


*0,a  — 


F  •  ka  •  exp 


(1  -  P)  ■  Ma  •  F  ■  E° 

R  ■  T 


*o,c  =  nc  ■  F  ■  kc  •  exp 


-P-nc-  F  ■  E° 
R  ■  T 


In  ref.  [9],  an  exhaustive  model  of  the  polarization  curve  is 
described;  in  the  model  the  anodic  overvoltage  is  neglected,  and 
the  value  of  the  cathode  exchange  current  density  is  considered 
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z'o,c  =  0.01  A  cm-2;  it  looks  a  high  value,  and  it  could  be  due  to 
the  fact  that  the  cathode  is  fed  with  pure  oxygen  in  place  of  air. 

In  refs.  [10,11],  a  unique  value  is  considered:  z'o  = 
4.84  x  10  8  A  cm  2;  the  low  value  could  be  explained  with 
the  materials  of  the  electrodes,  different  from  the  usual  carbon 
cloth;  moreover,  in  the  papers  a  dependence  of  the  electrode 
activation  from  the  ionic  conductivity  of  the  membrane  elec¬ 
trolyte  is  described,  thus  introducing  a  further  functional  link. 
In  ref.  [12],  the  considered  exchange  current  density  value  is: 
/o  =  1  x  10-6  A  cm-2. 

Finally,  there  is  an  interesting  paper  of  Parthasarathy  et  al. 
[13],  where  the  modification  of  the  exchange  current  density 
value  with  the  cell  operation  temperature  is  described.  At  low 
current  density,  the  platinum  particles  in  the  catalyst  layer  are 
covered  by  oxidate  particles  (oxide-covered  conditions)  separat¬ 
ing  them  from  the  oxygen  reactant  particles;  therefore,  at  lower 
load  the  reaction  speed  decreases  and  thus  the  exchange  current 
density:  the  values  at  the  cathode  electrode  are  in  the  range  from 
z’o,c  =  6  x  10-5  A  cm-2  at  T-  30  °C  to  z'o,c  =  2.6  x  10-4  A  cm-2 
at V=70°C. 

2.2.  Cell  resistance 

The  ohmic  overvoltage  could  be  expressed  by  ?iohm  = 
r-z  =  (rel  +  rion)-/,  with  the  electronic  (rel)  and  ionic  (Pon)  con¬ 
tributes.  The  electronic  resistance  increases  with  the  operation 
temperature.  The  ionic  resistance  of  the  membrane  is  related  to 
the  operation  temperature  too,  but  especially  to  the  degree  of 
humidification  of  the  membrane:  in  fact,  the  ionic  conductiv¬ 
ity  of  the  Nafion  increases  with  the  membrane  humidification. 
Therefore,  the  ionic  conductivity  could  have  a  direct  relationship 
with  the  current  density  (the  load).  In  ref.  [6],  there  is  a  analytic 
polynomial  expression  of  the  cell  resistance  found  from  exper¬ 
imental  tests  made  with  a  single  PEM  cell  of  50.56  cm2  active 
area  and  Nafion  117  membrane: 

r  =  Yi  +  Y2  •  T  +  y3  •  I  (3) 

where  y\  =0.811488  £2  cm2,  y2  =  —  1.7696  x  10-3  £2  cm2  K-1 
and  y3  =4.0488  x  10-3  Q  cm2  A-1. 

As  we  see,  according  to  the  empirical  evaluation  the  cell 
resistance  has  to  decrease  with  the  operation  temperature,  and 
to  increase  with  the  current  density.  As  an  example,  the  value 
of  the  cell  resistance  r  with  T  =  343.15  K  (70  °C)  and  1=  10  A 
(i  =  0.2  A  cm-2)  is  in  the  order  of  r  =  0.245  £2  cm2. 

In  a  further  paper  [14],  Amphlett  et  al.  introduce  another 
analytical  expression  of  the  ionic  resistance,  deducted  from  the 
analysis  of  literature  data  related  to  different  cells  all  with  a 
Nafion  117  membrane: 

ion  _  181.6  •  [1  +  0.03  •  i  +  0.062  •  (7/303)2  •  ;2-5] 

Y  ~  (A.  -  0.634  -  3  •  0  •  exp[4.18  •  ((T  -  303)/r)]  ' tm 

(4) 

The  cell  resistance  is  related  to  the  membrane  thickness  tm ,  to  the 
operation  temperature,  to  the  current  density  and  to  the  degree 
of  humidification  of  the  membrane  X.  The  parameter  X  has  been 
introduced  by  Springer  et  al.  [9],  and  it  is  defined  as  the  ratio 


between  the  number  of  water  molecules  and  the  number  of  sul¬ 
fonic  groups  S03“  of  the  membrane. 

Considering  indicative  values  (T=70°C,  A=14  correspond¬ 
ing  to  the  situation  of  saturated  reactant  flows  [15])  the  ionic 
resistance  increases  with  the  load,  with  an  approximate  linear 
relationship  with  a  slope  8.29  x  10-4  £2  m2  A-1. 

In  ref.  [9],  an  experimental  procedure  devoted  to  the  deter¬ 
mination  of  the  degree  of  humidification  of  the  membrane  is 
described,  which  allows  the  evaluation  of  the  ionic  resistance  in 
the  case  of  Nafion  1 17:  it  decreases  with  the  membrane  humidifi¬ 
cation,  and  an  indicative  value  with  T-  80  °C  and  i  =  0.5  A  cm-2 
is  Pon  =  0.285  ^  cm2.  The  value  described  in  Amphlett  et  al. 
[6]  is  lower  because  it  is  evaluated  at  lower  current  density 
(0.2  A  cm-2)  and  temperature  (with  a  coherence  with  (3)). 

In  ref.  [5],  the  cell  resistance  is  evaluated  at  T=10°C  and 
i  =  1  A  cm-2,  with  a  value  of  r  =  0. 15  Q  cm2,  without  specifying 
the  Nafion  thickness  (we  suppose  Nafion  115;  otherwise,  in  case 
of  Nafion  117,  the  value  would  not  be  coherent  with  the  values 
reported  by  Amphlett  and  Springer,  because  it  would  mean  a 
lower  resistance  in  case  of  higher  current  density). 

In  ref.  [16],  considering  a  Nafion  115  membrane  at  differ¬ 
ent  pressures  and  temperatures,  in  case  of  p  =  1  atm  at  T- 50  °C 
the  resistance  is  evaluated  as  r  =  0.363  £2  cm2,  and  at  T-  70  °C 
as  r  =  0.238  £2  cm2;  so  also  in  this  analysis,  the  resistance  r 
decreases  with  temperature.  In  ref.  [17],  Nafion  112,  115  and 
117  membranes  are  analyzed  in  the  range  70-80  °C,  at  different 
pressures,  and  the  parameters  of  a  polarization  curve  similar  to 
the  one  proposed  by  Kim  et  al.  [16],  are  estimated;  in  particular, 
in  case  of  Nafion  115  operating  at  atmospheric  pressure,  at  70  °C 
the  resistance  is  evaluated  as  r  =  0.24  £2  cm2.  Pisani  et  al.  [18] 
derived  a  semi-empirical  equation  of  the  polarization  curve  take 
into  account  especially  the  cell  voltage  deterioration  deriving 
from  the  cathode  active  region;  they  used  the  equation  to  fit  the 
set  of  experimental  data  of  [17]:  concerning  the  cell  resistance, 
in  case  of  Nafion  115  fed  by  H2/air,  operating  at  atmospheric 
pressure,  at  70  °C,  r  =  0.24  Q  cm2. 

2.3.  Internal  current  density 

The  internal  current  density  is  referred  to  the  electrons 
transported  through  the  electrolyte  membrane  and  to  the  fuel 
crossover.  This  current  is  thus  active  even  in  open  circuit  con¬ 
ditions  (equilibrium  of  the  electrochemical  reactions),  and  is 
related  to  the  reduced  value  of  the  open  circuit  voltage  (at  the 
PEMFC  temperature  range,  in  the  order  of  0.9  V,  see  the  graphs 
below)  compared  to  the  reversible  voltage  (in  the  order  of  1 .2  V). 
Its  effect  is  introduced  in  term  of  an  increase  of  the  activation, 
ohmic  and  concentration  overvoltages,  and  it  is  well  known  that 
the  internal  current  in  has  to  be  reduced  in  order  to  increase  the 
cell  behavior.  In  the  PEMFC  literature,  this  parameter  is  usually 
neglected,  or  it  is  considered  as  a  fixed  value  as  zn  =  2  mA  cm-2 
[4] .  In  ref.  [  1 9] ,  in  a  paper  describing  the  performances  of  a  mem¬ 
brane  operating  at  temperatures  higher  than  100  °C,  an  internal 
current  density  value  of  in  =  1  mA  cm-2  is  reported.  The  same 
author,  in  2005  [20]  analyzed  another  membrane  (28  \xm  thick) 
and  the  internal  current  density  is  around  2  mA  cm-2 .  In  a  recent 
paper  [21],  the  authors  underline  the  sensitivity  of  the  model  of 
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the  polarization  curve  to  the  internal  current  density,  which  only 
affects  the  simulation  results  at  low  current  values. 

A  paper  where  it  seems  to  be  introduced  a  concept  related  to 
the  internal  current  density  has  been  written  by  Kim  et  al.  [16]. 
An  analytical  expression  of  the  polarization  curve  in  the  form 
V  =  £°  —  b-\og  i  —  r-i  —  m-exp(A-/)  is  reported:  the  third  term  has 
an  expression  causing  an  overpotential  at  open  circuit  condi¬ 
tions.  The  parameter  ra,  expressed  in  V,  seems  to  determine  an 
effect  similar  to  the  one  caused  by  the  internal  current  density 
in  (in  the  paper  it  is  specified  that  the  physical  meaning  of  the 
parameters  m  and  n  is  not  well  understood).  In  the  paper,  the 
parameter  m  is  associated  to  a  value  in  the  order  of  the  mV,  thus 
maybe  too  low  to  express  the  reduction  of  the  open  circuit  volt¬ 
age  due  to  the  internal  current  density  (in  the  order  of  0.15  V). 

The  dependence  of  the  internal  current  density  from  opera¬ 
tion  variables  such  as  the  cell  temperature  is  not  discussed  in 
any  paper.  A  functional  relation  seems  to  exists:  an  increase  of 
the  operation  temperature  determines  an  improvement  of  the 
activation  of  the  electrochemical  reactions,  and  thus  it  could  be 
expected  an  increase  of  the  internal  current  density  value. 

2.4.  Limiting  current  density 

The  limiting  current  density  is  a  parameter  linked  to  the  con¬ 
centration  overpotential  at  the  electrodes,  which  is  significant 
just  at  high  load  values.  In  the  literature,  this  parameter  is  not 
particularly  discussed.  In  ref.  [2],  a  typical  value  for  a  PEM  is 
reported:  /i  =  0.9Acm-2.  In  the  electrochemical  literature  [6], 
an  analytical  expression  is  described,  where  the  limiting  cur¬ 
rent  is  function  of  the  diffusion  coefficient  of  the  reactant  in  the 
electrode,  of  the  electrode  thickness  and  of  the  reactant  concen¬ 
tration  above  the  diffusion  in  the  electrode.  Pisani  et  al.  [18] 
fitted  the  experimental  data  of  [17]  and  in  case  of  Nation  115 
fed  by  H2/air,  operating  at  atmospheric  pressure,  at  70  °C  the 
limiting  current  density  is  evaluated  as  i\  =  0.810  A  cm-2. 

3.  The  analytical  expression  adopted  for  the  PEMFC 
polarization  curve 

The  polarization  curve  of  a  single  PEM  fuel  cell  can  be 
described  by  the  analytical  expression  (1): 

Vc  =  L  —  T]aci  Oohm  Oconc  (5) 

where  E  is  the  open  circuit  voltage,  ^act  the  activation  overpoten¬ 
tial  in  the  two  electrodes,  r) 0hm  the  ohmic  overpotential  (ionic 
and  electronic)  and  Oconc  is  the  concentration  overpotential  in 
the  two  electrodes. 

The  four  terms  on  the  right  side  of  (5)  are  discussed  below. 
3.1.  Open  circuit  voltage 


If  the  reactants  and  the  product  are  assumed  as  having  ideal 
gas  behavior  (valid  approximation  because  of  the  low  operating 
temperatures  and  pressures  [22]),  the  activity  reduces  to  the  ratio 
of  partial  pressures  [1,8],  and  the  expression  (6)  is  modified  in 
(7): 


-A  g(T)  RT  PH  2-P°o52 

- 1 - -  In - - 

2  •  F  2  •  F  /?h2o 


where  the  operation  parameters  are  the  temperature  and  the  par¬ 
tial  pressures  of  the  reactants  and  the  product  in  the  interface 
catalyst  layer-membrane.  The  procedure  adopted  to  evaluate 
the  partial  pressures  in  (7)  is  reported  in  Appendix  A. 


3.2.  Activation  overpotential 


The  analytical  relation  between  the  overpotential  and  the 
current  density  on  an  electrode  surface  is  expressed  by  the 
Butler- Volmer  Eq.  (8)  [1,8]: 


exp 


RT  ) 


If  written  in  explicit  form  with  the  electrode  overpotential, 
assuming  that  the  transfer  coefficients  on  an  electrode  are  equal 
[1,23],  the  expression  (8)  becomes: 


0  act,  el 
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The  total  activation  overpotential  in  the  cell  is  the  sum  of  the 
anode  and  the  cathode  contributions: 


RT  .  * 

?7act  —  7?act,c  T  7?act,a  =  ~  ~  ’  sinh 

cr  •  t 


RT 

H - -  sinh 

aa  •  F 


(10) 


The  highest  activation  overpotential  is  at  the  cathode,  due  to 
the  lower  value  of  the  exchange  current  density  z‘o,c  [1,6].  The 
anode  and  cathode  exchange  currents  are  functions  of  many 
variables:  materials  and  porosity  of  the  electrode;  concentra¬ 
tion,  distribution  and  dimension  of  catalyst  particles;  operating 
temperature.  The  variable  which  can  be  modified  during  oper¬ 
ation  is  the  temperature:  a  temperature  increase  has  a  positive 
effect  on  the  semi-reaction  activation  on  the  electrode  surface, 
increasing  therefore  the  exchange  current  density  and  reducing 
the  activation  overpotential  [1,8]. 


3.3.  Ohmic  overpotential 


The  ohmic  overpotential  has  two  contributions,  linked  to  the 
electronic  and  ionic  resistance  [1]: 


*7 ohm  =  r-i  =  rel  •  i  +  rion 


(11) 


The  open  circuit  voltage  is  usually  expressed  by  Eq.  (6)  [1 ,8] : 
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The  ionic  resistance  can  be  expressed  as  a  function  of  the  ionic 
conductivity  and  the  thickness  of  the  electrolyte  membrane: 


rion 


fin 

cr 


(12) 
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The  ionic  conductivity  is  a  function  of  many  variables:  cell  tem¬ 
perature,  degree  of  humidification  of  the  membrane  (function 
of  cell  temperature,  current  density,  reactants  temperature  and 
humidification).  A  semi-empirical  expression  has  been  proposed 
by  Amphlett  et  al.  [14],  already  reported  as  Eq.  (4).  Expression 
(4)  is  general,  linking  the  ionic  resistance  to  the  temperature,  the 
current  density,  the  membrane  thickness  and  the  water  content 
of  the  membrane  X. 

The  parameter  X  has  been  introduced  by  Springer  et  al.  [9], 
and  express  the  number  of  water  molecules  associated  to  a 
sulfonic  group  SO3-  in  the  membrane.  The  degree  of  humidifi¬ 
cation  of  the  membrane  is  linked  to  the  water  produced  by  the 
cathode  semi-reaction  (and  therefore  to  the  current  load),  to  the 
inlet  water  mixed  with  the  reactants  (that  is,  to  the  degree  of 
humidification  of  the  reactant  gases),  and  to  the  mass  transport 
phenomena  occurring  in  the  membrane  [9,15,24-27]. 

Other  models  are  reported  in  literature  to  describe  the  cell 
resistance  [6,14]. 


3.4.  Concentration  overpotential 


The  concentration  overpotential  on  a  electrode  can  be 
expressed  as  [1,8,10,28]: 
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often  expressed  introducing  the  concept  of  limiting  current  den¬ 
sity: 
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The  total  concentration  overpotential  of  the  cell  is  the  sum  of 
the  anode  and  cathode  contributions: 
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3.5.  Internal  current 

Even  at  open  circuit,  the  experimental  observation  shows  that 
the  reversible  voltage  (7)  is  not  reached  by  the  cell.  At  open  cir¬ 
cuit,  the  anodic  semi-reaction  of  hydrogen  is  at  equilibrium, 
but  it  could  happen  that  a  very  limited  fraction  of  the  hydro¬ 
gen  molecules  oxidize  producing  H+  ions  and  electrons;  both 
diffuse  through  the  membrane  (there  is  not  an  external  load  con¬ 
nected)  to  the  cathode  side,  generating  a  current.  At  the  same 
time,  there  is  the  possibility  that  a  very  limited  part  of  the  hydro¬ 
gen  molecules  does  not  participate  to  the  anodic  semi-reaction, 
and  diffuse  to  the  cathode  in  form  of  molecules  (fuel  crossover). 
The  two  phenomena  are  modeled  as  a  current  flow  inside  the 
membrane,  introducing  the  concept  of  internal  current  in.  The 
internal  current  could  occur  also  at  open  circuit  conditions,  and 
therefore  it  could  explain  the  voltage  reduction  compared  to  the 


reversible  voltage.  Therefore,  every  value  of  the  current  density 
in  Eqs.  (10),  (11)  and  (15)  is  increased  by  the  internal  current 
term  in  [4] . 

3.6.  Analytical  expression  of  the  PEMFC  polarization 
curve 


The  complete  analytical  expression  of  the  PEMFC  polariza¬ 
tion  curve  is: 
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(16) 


The  expression  (16)  will  be  used  to  fit  the  experimental  polar¬ 
ization  curves  described  in  the  following. 

4.  Parameters  estimation  of  the  polarization  curve 


In  the  analytical  expression  of  the  PEMFC  polarization  curve, 
Eq.  (16),  the  independent  variables  are  the  current  density,  the 
temperature  and  the  pressure,  while  the  dependent  variable  is 
the  voltage.  Moreover  there  are  a  few  constants,  that  usually 
are  called  parameters,  and  some  of  these  are  often  unknown. 
When  an  experiment  is  done  to  obtain  the  polarization  curve, 
the  voltage  (dependent  variable)  is  measured  directly,  but  the 
parameters  are  not. 

The  estimation  of  some  of  these  parameters  is  an  aim  of  this 
work.  The  problem  is  to  obtain  the  good  estimation  of  these 
parameters  using  a  discrete  number  of  experimental  results,  but 
often  not  all  the  parameters  can  be  independently  estimate.  It  is 
useful  to  rewrite  Eq.  (16)  in  the  form: 

V  =  /(/,  T,  p ,  aa,  ac,  z'o, a,  *o,c  h, a,  h,c  r,  in)  (17) 

In  Eq.  (17),  eight  parameters  are  pointed  out.  For  some  of  these 
parameters  further  consideration  have  been  done. 

Concerning  the  transfer  coefficient  on  electrode,  it  can  be 
defined  as  [8]: 


=  (1  -  P)  -  ftel 
ac  =  P  •  nQ  1 


(18) 


where  nQ\  is  the  number  of  electrons  in  the  rate  determining  step 
of  the  reaction  (4  for  the  anode  and  1  for  the  cathode)  and  ft  is  a 
symmetry  factor  whose  value  is  very  near  to  0.5  [1 ,23] ;  therefore, 
for  the  transfer  coefficients  the  values  aa  =  2  and  ac  =  0.5  have 
been  assumed. 

Concerning  the  limiting  current  densities,  their  values  are 
estimated  to  be  very  high  (ca.  43  A  cm-2  for  anode  and  9  A  cm-2 
for  cathode),  but  in  the  cell  analyzed  in  this  paper  the  maximum 
current  density  imposed  is  0.25  A  cm-2,  because  the  adopted 
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ME  A  (ElectroChem  Inc.)  showed  a  sudden  voltage  drop  due  to 
the  low  ionic  conductivity  of  the  membrane.  So,  we  conclude 
that  the  effect  of  the  concentration  losses  is  negligible  and  in  our 
analysis  the  parameters  z^a  and  i\?c  are  not  meaningful. 

Concerning  the  area  specific  resistance  r,  it  could  be  divided 
in  two  components 

r  =  rel  +  rion  (19) 

where  rel  is  the  electronic  component  and  Aon  is  the  ionic.  The 
rel  is  function  of  the  temperature,  while  for  Pon  Amphlett  et  al. 
[14]  proposed  the  general  law  [14,29]  already  reported  in  Eq. 
(4).  Fixing  X  and  tm,  Eq.  (4)  could  be  approximated  as: 

rion  =  C0(T)  +  C2-i  (20) 

therefore 

r  =  rtX(T)  +  C0(T)  +  C2  ■  i  =  CX{T)  +  C2  ■  i  (21) 

Summarizing,  Eq.  (17)  can  be  now  rewritten  in  the  form: 

V  =  /(/,/ofa,io,c,C1,C2,in)  (22) 

with  five  parameters  unknown. 

To  verify  which  parameters  can  be  simultaneously  estimated 
it  is  useful  to  analyze  the  sensitivity  coefficients:  the  ith  sen¬ 
sitivity  coefficient  of  a  parameter  is  the  first  derivate  of  the 
function  (22)  with  respect  to  the  ith  parameter,  multiplied  by 
the  ith  parameter  itself  [3]: 

dv 

Xi  =  —Pi  (23) 

0Pi 

When  these  coefficients,  over  the  range  of  the  observation,  are 
not  linearly  dependent  all  parameters  can  be  estimated  simulta¬ 
neously. 

To  verify  if  the  linear  dependence  occurs,  these  coefficients 
were  plotted  (Fig.  1):  there  is  a  linear  dependence  between  C\ 
and  C2  and  between  C\  (or  C2)  and  z'o,a-  Therefore,  we  decided 
to  estimate  the  parameter  z'o,c,  C\  and  in.  For  the  analyzed  tem¬ 
perature  range,  the  values  of  C2 ,  evaluated  through  Eq.  (4),  are 
shown  in  Table  1,  where  from  Weber  and  Newman  [15]  a  value 
of  X  =  14  has  been  assumed,  valid  for  a  membrane  in  contact 
with  saturated  water,  and  for  the  used  membrane  the  thickness 
was  127  (Jim. 
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Fig.  1.  Sensitivity  coefficient  vs.  current  density. 


Table  1 


Values  of  the  parameter  C2  (X  =  14,  membrane  thickness  127  |xm  in  contact  with 
saturated  water)  in  the  temperature  range  50-80  °C 


r(°c) 

C2  (92  cm4  A”* 1) 

50 

0.0416 

55 

0.0393 

60 

0.0372 

65 

0.0353 

70 

0.0335 

75 

0.0318 

80 

0.0303 

%  error  on  i0a 

Fig.  2.  Relative  error  on  estimated  parameters  due  to  an  imposed  error  on  *o,a- 

The  exchange  current  density  at  the  anode  z'o,a  was  assumed 
from  literature  [4]  equal  to  0.2  A  cm-2.  To  verify  the  effect  on 
the  parameter  estimation  of  a  deviation  from  this  value,  an  anal¬ 
ysis  using  experimental  simulations  was  done.  That  is,  first  some 
polarization  curves  were  simulated  using  different  value  of  z'o,a, 
then  the  other  parameters  were  estimated  imposing  always  the 
reference  value  of  z'o,a.  An  analysis  imposing  a  known  error 
on  /o, a  (±25%)  and  considering  the  two  extreme  temperatures 
(50  °C  and  80  °C)  was  done.  The  results  of  this  analysis  are 
shown  in  Fig.  2:  the  greater  influence  was  for  the  parameter  zb,c  at 
low  temperature,  in  fact  the  maximum  deviation  was  10.4%.  At 
high  temperature  the  maximum  effect  was  less  important  (3.1%). 
Also  the  parameter  in  was  influenced  from  an  zb, a  error  especially 
at  lower  temperature  (4.0%,  compared  to  2.5%  at  higher  tem¬ 
perature).  On  parameter  C\  there  was  a  lower  influence  when 
the  temperature  was  high  (1.3%)  while  if  low  temperature  was 
analyzed  the  maximum  error  was  very  little  (0.3%). 

5.  Experimental 

The  fuel  cell  employed  for  the  experiments  is  an  Elec¬ 
troChem  Inc.  single  cell  (EFC25-01SP).  This  cell  employs  a 
Nation  115  membrane  (127  pan  thickness),  with  an  active  area 
of  25  cm2 *.  The  electrodes  are  identical  and  use  Toray™  ca. 
0.20  mm  thick  carbon  paper  for  the  backing;  the  catalyst  layer  is 
on  the  order  of  ca.  0.05  mm  thick  so  the  total  electrode  thickness 
is  0.25  ±  0.02  mm.  The  catalyst  layer  is  made  of  carbon  sup¬ 
ported  platinum,  loading  lmgcm-2  20wt.%  Pt/C.  The  mem¬ 
brane  electrode  assembly  (MEA)  was  first  sandwiched  between 
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|  Water  make-up  Anode  Inlet 


Fig.  3.  Schematic  drawing  of  the  gas  flow  apparatus. 


two  gaskets  and  then  between  two  graphite  current  collector 
plates. 

The  analyzed  fuel  cell  has  been  tested  with  a  station  made 
by  ElectroChem  Inc.  (named  Power  Station  CompuCell );  with 
this  test  station  three  different  fluids  was  managed:  the  oxidant 
gas  (in  this  study  air),  the  fuel  (hydrogen)  and  an  inert  gas  for 
purging  purposes  (nitrogen).  The  gas  flow  rates  are  controlled 
with  mass-flow  controllers.  The  Gas  Management  Unit  (GMU) 
can  supply  humidified  or  dry  gas  to  the  fuel  cell.  The  humidity 
is  imposed  by  sparging  the  gases  through  two  gurgle  tanks  filled 
with  distilled  water.  The  tanks  are  heated  and  kept  at  the  chosen 
temperatures  by  controllers.  The  temperature  sensors  are  K-type 
thermocouples.  The  scheme  of  the  experimental  apparatus  is 
shown  in  Fig.  3. 

A  set  of  experimental  tests  has  been  scheduled  to  analyze  the 
effects  of  the  temperature  on  cell  performance  and  especially 
on  the  behavior  of  some  characteristic  parameters  (exchange 
and  internal  current  densities;  membrane  resistance).  To  out¬ 
line  the  temperature  effect,  the  tests  have  been  carried  out 
imposing  the  cell  temperature  in  the  range  T  =  50-80  °C  (max¬ 
imum  design  temperature  of  the  experimental  setup),  with  a 
temperature  step  of  5  °C  (with  7  runs  for  every  session);  the 
number  of  sessions  has  been  5,  and  therefore  the  total  number 
of  runs  has  been  35.  The  same  temperature  has  been  adopted 
for  both  the  inlet  streams,  and  moreover  they  were  in  satura¬ 
tion  conditions,  therefore  T  =  Tas  =  Tcs.  The  tests  have  been 
performed  at  ambient  pressure  (reactant  backpressures  fixed  to 
zero). 


The  available  test  facility  does  not  allows  an  automatic  mod¬ 
ification  of  the  reactant  mass  flows  according  to  the  load,  and 
therefore  operates  with  a  constant  mass  flow  during  the  runs. 
The  maximum  value  of  the  cell  current  density  has  been  fixed  at 
0.25  A  cm-2  (7.5  A),  but  the  volume  flow  values  were  fixed  in 
order  to  reach  a  higher  current  load  (0.8  A  cm-2):  141  ml  min-1 
of  hydrogen  (excess  Xu  =  1.2)  and  726  ml  min-1  of  air  (excess 
A. air  =  2).  As  these  volume  flows  are  fixed,  a  further  excess  of 
hydrogen  and  air  at  low  current  densities  occurs. 

Some  preliminary  experiments  were  carried  out  to  choose  the 
load  profile  that  was  then  adopted  in  the  tests.  As  observed  above, 
the  maximum  current  density  was  fixed  at  0.25  A  cm-2 .  The  cur¬ 
rent  density  therefore  varies  between  0  A  cm-2  and  0.25  A  cm-2 
with  steps  of  0.025  A  cm-2 .  Every  current  density  plateau  is  kept 
for  30  s,  and  the  ramp  between  the  plateau  also  lasts  30  s.  There¬ 
fore,  the  total  test  length  was  ca.  1 1  min.  The  sampling  rate  for 
data  logging  was  1  s.  The  considered  voltage  was  the  average 
value,  and  the  estimated  maximum  voltage  uncertainty  is  less 
than  3%. 

6.  Results  and  comments 

6.1.  Regression  analysis 

Fig.  4  shows  how  the  model  can  rebuild  the  polarization  curve 
of  the  cell  when  the  unknown  parameters  are  estimated  at  a 
fixed  temperature.  In  the  figure  the  regression  curves  obtained 
at  different  values  of  the  temperature  range  (lower  bound: 
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Fig.  4.  Regression  curves  obtained  at  different  values  of  the  temperature 
range — 50  °C:  rhombus  experimental,  solid  line  model;  65  °C:  asterisks  exper¬ 
imental,  dashed  line  model;  80  °C:  circles  experimental,  dotted  line  model. 
Maximum  voltage  uncertainty  less  than  3%. 

50  °C;  intermediate  value:  65  °C;  upper  bound:  80  °C)  are 
shown. 

As  it  was  expected,  the  activation  of  the  reactions  is  facilitated 
at  higher  temperatures  (lower  activation  overpotential),  and  the 
cell  resistance  is  lower  at  higher  temperatures  (lower  slope  of 
the  curve).  From  the  graph  it  is  possible  to  see  that  the  analytic 
curves  rebuild  the  experimental  markers  very  satisfactorily.  We 
only  notice  a  bigger  distance  between  the  markers  and  the  curves 
at  0.025  A  cm-2 :  we  think  this  fact  is  due  to  the  estimation  of  the 
cathode  exchange  current  density,  which  is  the  parameter  most 
sensitive  to  errors  and  the  most  important  parameter  for  the  part 
of  the  curve  where  the  imprecision  has  been  noticed. 

Concerning  the  performance  of  the  adopted  MEA  (Elec- 
troChem  Inc.),  it  is  evident  that  the  cell  overpotentials  are  very 
high:  at  a  low  value  of  current  density  (e.g.  0.2  A  cm-2)  the  cell 
voltage  is  already  dropped  at  values  around  0.4  V.  As  shown  in 
Fig.  4,  the  sudden  voltage  drop  is  mainly  due  to  a  low  ionic 
conductivity  of  the  membrane,  and  this  effect  is  confirmed  by 
the  high  values  obtained  by  the  parameter  C\  of  Eqs.  (20)  and 
(21)  modeling  the  cell  resistance,  which  are  discussed  below. 
Nevertheless,  the  analysis  of  the  performances  of  the  MEA  do 
not  represent  the  aim  of  the  paper,  which  concentrates  on  the 
discussion  of  the  MEA  model  and  of  the  regression  procedures 
applied  to  experimental  data  to  estimate  some  important  cell 
parameters. 

6.2.  Dependence  ofio,c  with  temperature 

In  Fig.  5,  the  estimated  values  and  the  error  bars  of  the  param¬ 
eter  z’o,c  versus  the  operating  temperatures  are  shown. 

Concerning  the  obtained  values,  they  are  coherent  with  other 
values  indicated  in  literature  [4-13].  In  particular,  the  trend  of 
the  parameter  with  temperature  seems  very  similar  to  the  trend 
reported  by  Parthasarathy  et  al.  [13],  and  this  could  confirm  the 
validity  of  the  adopted  model  and  especially  of  the  parameter 
estimation  procedure.  From  the  graph,  it  is  possible  to  observe  a 
qualitative  behavior  of  the  parameter:  it  increases  with  tempera¬ 
ture,  and  the  slope  of  the  function  is  close  to  a  constant  value  in 


xIO-4 


the  range  50-75  °C  (nearly  a  linear  shape);  then  a  slope  incre¬ 
ment  occurs  at  greater  temperatures.  The  absolute  values  of  the 
parameter  are  lower  than  the  ones  obtained  by  Parthasarathy  et 
al.  [  1 3] ,  and  this  is  linked  to  the  observed  low  quality  of  the  MEA; 
nevertheless,  the  values  are  of  the  same  order  of  magnitude 
(10-4  A  cm-2),  marking  that  the  cathode  activation  overpoten¬ 
tial  does  not  represent  the  main  cause  of  the  low  performance  of 
the  MEA.  The  estimated  uncertainty  of  this  parameter  is  ±10%. 

6.3.  Dependence  of  r  (Cj)  with  temperature 

In  Fig.  6,  the  estimated  values  and  the  error  bars  of  the  param¬ 
eter  Ci  versus  the  operating  temperatures  are  shown. 

The  trend  of  the  parameter  with  temperature  is  the  expected 
one,  typical  of  the  electrolytes  (ionic  conductors):  the  qualitative 
behavior  shows  a  quite  linear  decrease  in  the  range  of  tempera¬ 
tures  from  50  °C  to  80  °C.  The  absolute  values  of  the  parameter 
are  consistently  higher  than  the  ones  reported  in  literature  related 
to  other  ME  As:  the  obtained  value  is  in  the  range  1 .6-1 .9  ^  cm2 
(estimated  uncertainty  is  less  than  ±1%),  while  typical  values 
are  in  the  order  of  0.3  £2  cm2  [5,6,9,14,16].  As  discussed  above 
in  case  of  z'o,c>  this  is  linked  to  the  evident  low  quality  of  the 
tested  MEA:  in  particular,  these  values  of  cell  resistance,  and 
the  related  high  values  of  ohmic  overpotentials,  seem  to  repre¬ 
sent  the  main  cause  of  the  low  performance  of  the  MEA. 


Fig.  6.  Values  of  the  parameter  C\  vs.  cell  operating  temperature. 
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Fig.  7.  Values  of  the  parameter  in  vs.  cell  operating  temperature. 


6.4.  Dependence  ofin  with  temperature 

In  Fig.  7,  the  estimated  values  and  the  error  bars  of  the  param¬ 
eter  in  versus  the  operating  temperatures  are  shown. 

As  stated  above,  in  literature  we  have  not  found  descrip¬ 
tions  of  a  possible  trend  of  this  parameter  with  temperature.  The 
trend  observed  with  our  data  shows  an  increase  of  the  param¬ 
eter  with  temperature  in  the  range  from  50  °C  to  80  °C.  The 
absolute  values  of  the  parameter  are  higher  than  the  only  case 
reported  in  literature  [4]:  the  obtained  value  is  in  the  range 
6-12  mA  cm-2  (estimated  uncertainty  is  ca.  ±7%),  compared 
with  the  value  indicated  by  Larminie  and  Dicks  which  is  in 
the  order  of  2mA  cm-2.  Apparently,  these  values  seem  to  rep¬ 
resent  another  cause  of  the  low  performance  of  the  adopted 
MEA,  because  a  high  value  of  in  means  an  increase  of  the  fuel 
crossover.  At  the  same  time,  the  trend  increasing  with  tempera¬ 
ture  seems  to  suggest  that  the  crossover  phenomena  occurs  with 
higher  probability  at  increased  temperatures.  The  last  sentences 
have  to  be  analyzed  in  more  detail  with  other  experimental  ses¬ 
sions. 

7.  Application  of  the  regression  model  to  the 
experimental  data  of  other  authors  and  discussion 

To  give  some  evidence  of  the  validity  of  the  regression  model 
and  procedure,  they  have  been  extended  to  other  sets  of  experi¬ 
mental  data,  besides  the  experimental  sessions  produced  in  our 
laboratory.  Experimental  data  in  form  of  polarization  curves  of 
single  cells  are  largely  available  in  literature.  Unfortunately,  not 
many  papers  report  also  the  results  of  the  fitting  procedure  in 
terms  of  the  values  of  the  parameters  of  the  polarization  curve 
with  their  uncertainty;  these  values  are  necessary  in  order  to 


check  the  regression  model  and  procedure.  We  have  considered 
in  particular  two  papers  [16,17].  Both  the  papers  report  experi¬ 
mental  data  in  form  of  polarization  curves,  and  especially  report 
the  estimated  values  of  one  of  the  parameters  considered  in  our 
paper:  the  cell  resistance  r  (actually,  our  model  has  separated  the 
cell  resistance  in  an  expression  with  two  parameters,  Eq.  (21), 
but  from  the  values  of  C\  and  C2  we  can  obtain  r). 

As  shown  in  Table  2,  the  values  of  r  found  with  our  regression 
model  are  in  good  agreement  with  the  values  of  the  parameter 
reported  by  Kim  et  al.  The  values  estimated  by  our  regression 
model  are  slightly  lower  than  the  values  reported  by  Kim  et  al.: 
therefore,  we  assign  a  lower  value  to  the  ohmic  overpotential.  A 
possible  explanation  can  be  linked  to  the  role  and  values  of  the 
cathode  exchange  current  density  parameter  z'o,c-  We  have  tried 
to  estimate  the  values  of  the  parameter  from  the  data  available 
by  the  paper.  The  model  indicates  as  Eo  =  Ex  +  Mog  io,  with  Eo 
is  the  fitted  parameter,  Er  the  reversible  potential  for  the  cell, 
and  io  and  h  are  the  Tafel  parameters  for  oxygen  reduction.  The 
values  of  Eo  and  h  are  reported  by  Kim  et  al.;  the  value  of  Ev 
can  be  read  in  Figs.  4  and  5  of  [16]  (with  some  approxima¬ 
tion);  finally,  the  value  of  io  can  be  estimated  in  the  order  of 
5-6  x  10-1  A  cm-2.  These  values  are  higher  than  other  values 
reported  in  literature  (in  the  order  of  1-2  x  10  4  A  cm  2)  and 
even  than  the  values  estimated  by  our  regression  model  (in  the 
order  of  3-4  x  10-4  A  cm-2).  This  means  that,  according  to  the 
estimated  values,  the  cathode  activation  overpotential  are  very 
low  in  the  regression  model  used  by  Kim  et  al.  This  causes  the 
fact  that  their  regression  model  has  to  assign  a  greater  weight 
to  the  ohmic  overpotential,  explaining  the  higher  values  of  the 
parameter  r  compared  by  the  values  found  by  our  regression 
model.  Maybe  the  regression  model  used  by  Kim  et  al.  underes¬ 
timates  the  weight  of  the  cathode  activation  overpotential,  which 
is  not  negligible  compared  to  the  ohmic  overpotential  (see,  for 
example  [31],  Fig.  4,  p.  2483). 

Concerning  the  behavior  of  r  versus  T  the  slope  is 
5.7  x  10-3  Q  cm2  K-1  for  our  experimental  results  and  a  slope 
of  6.2  x  10  3  £2cm2K  1  for  Kim  et  al.  experimental  results. 
These  results  are  in  good  agreement. 

Concerning  the  paper  of  Squadrito  et  al.  [17],  the  discus¬ 
sion  is  similar.  As  shown  in  Table  2,  the  values  of  r  found 
with  our  regression  model  are  lower  than  the  values  reported 
by  Squadrito  et  al.  But,  as  reported  in  Squadrito  et  al.  (p.  1454), 
the  expected  values  were  in  the  order  of  0.12-0.14  Q  cm2,  in 
accordance  with  the  value  estimated  by  our  regression  model 
applied  to  the  data  provided  by  the  authors.  In  our  opin¬ 
ion,  even  in  this  case  a  possible  explanation  can  be  linked  to 
the  role  and  values  of  the  cathode  exchange  current  density 
parameter. 


Table  2 

Comparison  between  the  values  of  the  cell  resistance  evaluated  with  our  regression  model  and  the  values  reported  by  Kim  et  al.  [16]  and  Squadrito  et  al.  [17] 


Cell  temperature  (°C) 

Data  from  Kim  et  al.  [16] 

Data  from  Squadrito  et  al.  [17] 

Eq.  (16),  our  paper  (£2  cm2) 

Eq.  (5)  [16]  (12  cm2) 

Eq.  (16),  our  paper  (£2  cm2) 

Eq.  (5)  [17]  (12  cm2) 

50 

0.323 

0.363 

— 

— 

70 

0.208 

0.238 

0.135 

0.240 

Membrane:  Nation  115;  reactant:  tT/air;  pressure:  1  atm. 
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8.  Conclusions 


catalyst  layer)  the  Stefan-Maxwell  equation  has  been  used  [34] : 


The  paper  describes  a  methodology  for  the  estimation  of  some 
parameters  of  the  model  of  the  polarization  curve  of  the  cell,  and 
analyzes  and  discusses  the  behavior  of  the  parameters  versus  a 
fundamental  independent  variable  (the  cell  operating  tempera¬ 
ture).  The  results  can  be  summarized  as  follows: 

•  from  the  parameters  analysis  it  has  been  verified  that  only 
three  parameters  of  the  cell  polarization  curve  model  can  be 
simultaneously  estimated:  /o,c>  C\  and  /n;  the  low  sensibility 
of  the  model  due  to  the  effects  of  the  other  parameters  has 
been  demonstrated  using  experimental  simulations; 

•  the  analytical  model  of  the  polarization  curve  rebuilds  the 
experimental  results  very  satisfactorily;  the  regression  is  less 
accurate  at  low  current  densities,  probably  due  to  the  esti¬ 
mation  of  the  cathode  exchange  current  density  which  is  the 
parameter  most  sensitive  to  errors; 

•  the  qualitative  behavior  of  the  parameter  z'o,c  shows  a  quite 
linear  increase  in  the  range  of  temperatures  from  50  °C  to 
80  °C,  with  values  in  coherence  with  the  literature  reports 
(10-4  A  cm-2); 

•  the  trend  of  the  parameter  C\  (linked  to  the  cell  resistance)  is 
typical  of  the  electrolytes,  with  a  quite  linear  decrease  with 
temperature; 

•  the  qualitative  behavior  of  the  parameter  in  shows  an  impor¬ 
tant  increase  in  the  range  of  temperatures  from  50  °C  to 
80  °C;  the  increase  with  temperature  seems  to  suggest  that 
the  crossover  phenomena  occurs  with  higher  probability  at 
increased  temperatures; 

•  due  to  the  limited  numbers  of  experimental  points,  and  espe¬ 
cially  to  the  fact  that  the  experiments  have  been  done  using 
just  one  ME  A,  it  is  not  possible  to  propose  general  analytic 
correlations  linking  the  parameters  and  the  operating  tem¬ 
perature,  but  the  potentiality  of  the  method  of  analysis  and 
the  easiness  of  the  experiments  allow  to  apply  it  for  a  wide 
typologies  of  ME  As;  in  fact  we  have  checked  the  regression 
model  with  data  of  cell  resistance  provided  by  other  authors 
with  good  agreement  of  the  results. 
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where  N  is  the  molar  flux  of  a  component  per  unit  surface. 
The  binary  diffusion  coefficient  D/y  have  been  evaluated  through 
expression  (A2): 
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where  the  coefficients  a  and  b  depend  on  the  mixture  and  are 
reported  in  [34]. 

The  effective  diffusivity  in  the  porous  diffusion  layer  has  to 
be  evaluated  from  the  binary  diffusion  coefficient.  There  are 
different  models  proposed  to  evaluate  the  effect  of  porosity  s 
of  the  transport  medium  [5].  We  have  considered  the  correla¬ 
tion  suggested  by  [35],  valid  for  random  fibrous  porous  media 
representing  the  diffusion  layer  of  the  cell  electrode: 


Deff  =  £>a,b  •  s  • 


Y 


(A3) 


Concerning  the  porosity,  in  absence  of  specific  data  for  our  cell, 
we  have  considered  a  value  widely  used  in  literature:  s  =  0.3  [6]. 
The  parameter  £p  is  a  percolation  threshold:  for  porous  media 
composed  of  two-dimensional,  long  and  overlapping  random 
fibber  layers,  it  has  been  found  as  0.11  [5].  The  term  y  is  an 
empirical  constant,  found  as  0.785  in  case  of  cross-plane  diffu¬ 
sion  [5]. 

The  molar  flux  per  surface  area  N  represents  the  net  number 
of  moles  of  a  chemical  species  transported  through  the  diffu¬ 
sion  layer  to  the  membrane  interface  of  active  surface  S  for  unit 
time.  The  number  of  moles  participating  to  the  reaction  for  unit 
time  and  surface  is  evaluated  with  the  Faraday  law — concerning 
hydrogen  and  oxygen  we  have: 
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Appendix  A.  Procedure  adopted  to  evaluate  the  partial 
pressures  in  Eq.  (7) 


Nitrogen  is  inert,  and  thus  its  net  flux  through  the  diffusion  layer 
is  An2  =  0.  Concerning  the  water  molar  flux,  we  have  to  con¬ 
sider  both  the  inlet  water  (with  humidified  reactant  flows)  and 
the  transport  mechanisms  inside  the  membrane.  The  electro- 
osmotic  drag  of  water  from  anode  to  cathode  is  described 
through  [24,36]: 


The  partial  pressures  at  the  interface  are  evaluated  through 
laws  of  diffusion  in  porous  media  [5,29-33].  Considering  the 
flow  in  the  inlet  ducts  of  the  bipolar  plate,  a  mass  fraction  diffuses 
through  the  electrode  diffusion  layer.  The  mass  transport  has 
been  considered  perpendicular  to  the  diffusion  layer  surface.  To 
evaluate  the  molar  fraction  y  of  each  component  of  the  mixed 
flow  after  the  diffusion  layer  (that  is,  in  correspondence  of  the 


(A5) 


where  n&  is  the  electro-osmotic  drag  coefficient.  In  literature, 
different  values  of  n&  can  be  found,  directly  linked  to  the  degree 
of  humidification  of  the  membrane  [9,15,36];  using  the  exper¬ 
imental  data  described  in  [24],  we  have  estimated  the  value 
n&  =  0.27  molH2o  mol^j  •  The  other  transport  mechanism  is  the 
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bask  diffusion  from  cathode  to  anode,  linked  to  the  concentra¬ 
tion  gradient  at  the  two  sides  of  the  membrane  (with  higher  value 
at  the  cathode).  The  molar  flux  of  water  due  to  back  diffusion  is 
expressed  by  [9,34,36]: 


^H20,bd  =  Av 


CH2Q,c  <Yl20,a 

6n 


(A6) 


The  water  diffusion  coefficient  Dw  is  evaluated  through  expres¬ 
sions  found  in  [9]  as  Dw  =  1.28  x  10“ 10  m2  s-1. 

The  molar  flux  of  water  to  the  anode  reaction  site  coming 
from  the  channels  (from  the  humidified  reactant  flow)  is  Aht2o,a; 
a  fraction  migrates  to  the  cathode  via  electro-osmosis  drag,  while 
another  molar  flux  diffuses  to  the  anode  via  back  diffusion  from 
the  cathode  reaction  site;  the  net  water  molar  flux  through  the 
membrane  (positive  if  directed  to  cathode)  is  therefore: 


AH20,net  =  An20,eo  —  An20,bd  (A7) 

The  water  molar  flux  which  returns  from  the  anode  reaction  site 
to  the  anode  channels  is  thus  A^H2o,a  —  AH2o,net-  So,  the  net 
water  molar  flux  transported  through  the  anode  diffusion  layer 
is  evaluated  as: 


AH20,a,net  =  An20,a  “  (An20,a  —  An20,net)  =  An20,net 


(A8) 

At  the  cathode,  Aht2o,c  is  the  molar  flux  of  water  to  the  cathode 
reaction  site  from  the  channels  (from  the  humidified  reactant 
flow);  A^H2o,net  is  the  molar  flux  coming  from  the  membrane; 
moreover,  A^H2o,prod  is  the  molar  flux  of  produced  water.  There¬ 
fore,  the  net  water  molar  flux  transported  through  the  cathode 
diffusion  layer  is: 

AH20,c,net  “  (Ar20,c  +  AH20,prod  “  An20,net) 

—  —  (An20,net  —  An20,prod)  (A9) 

where  the  negative  sign  stands  for  a  flux  from  the  catalyst  layer 
to  the  distribution  channels.  There  are  all  the  details  to  solve  the 
differential  equation  (A2). 
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